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1,  INTRODUCTION 


^  /->v* 


In  this  report  we  give  finite  difference  approximations 
to  the  differential  equations  for  the  following  problems: 

(1)  the  bending  of  a  variable -thickness  truncated  right  coni¬ 
cal  shell  attached  to  an  axially  and  sectlonally  symmetric 
base  plate  of  variable  thickness,  (2)  the  bending  of  two 
Joined  variable -thickness  truncated  right  conical  shells. 

In  case  (l)  above,  the  terms  arising  from  consideration 
of  the  transverse  deflection  of  the  base  caused  by  radial  load¬ 
ing  have  been  included.  As  a  consequence  the  finite  differ¬ 
ence  equations  are  non-linear.  For  solutions  of  these  equa¬ 
tions  we  suggest  an  iterative  s±eme  combined  with  a  direct 
method  based  on  the  factorization  theorem  for  square  matrices. 
In  case  (2),  the  finite  difference  equations  are  linear  and 
procedures  for  solution  are  straightforward. 

II .  DIPPEHEf'ITIAL  EQUATIONS 

The  differential  equations  for  the  conical  shells  are 
derived  in  TR  1-21  [1]  and  are  repeated  below  for  convenience. 
Coordinate  system,  geometry,  and  variables  are  shox«i  in  Pig.  1. 
The  differential  equations  are  given  in  terms  of  the  variables 
S  and  6  where  S  is  related  to  the  force  per  unit  length 
(see  Pig.  2)  by: 

Numbers  in  square  brackets  refer  to  bibliography  on  Page  29« 
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2.1 


and  8  Is  the  slope  of  the  median  plane  in  the  axial  direction, 
l.e. : 


The  differential  equations  are  ([11  p.  l6): 

2.3  L(g)  +  f^9  =  +  F(s) 

2.4  L(S)  +  fgS  ==  XgB 

where  the  differential  operator  L  Is  of  second  order: 


L(ll)  =  h  cot^fau"  +  (1  +  3s  |-)u'  -  I" 


The  apostrophe  rex)resents  differentiation  with  respect  to  s. 
The  other  quantities  appearing  in  the  differential  equations 
are  defined  as  follows: 


f^  =  svh  cot*^ 


where  v  is  Poisson's  ratio. 


«  4  - 


2.7 


fg  *=  f  (2  +  v)h*  +  2sh  ^  cot<^ 

2.8  = 

where  E  is  the  modulus  of  elasticity, 

2.9 
and 

2.10 


Xg  —  ~  K 

F(s)  = 

^  h*^  cotcp 


where  ([11  p.  4): 


2.11 


s 

p(8)  =  -  fsN  slnf,5  +  s  Q  coscp  !  =  {  Z  cospp  sds 
-  *0  J  J« 

*^0 


m  sin?p  +  sQ.  cos'fl 


'J, 


S”S, 


and  Z  is  a  distributed  load  normal  to  the  middle  surface  of 
the  shell  and  acting  over  the  whole  surface,  and  is  a  shear 
force  per  unit  length  acting  normal  to  the  middle  surface  of 
the  shell  (see  Fig. 2) 

The  differential  equations  for  the  base  may  be  written 


-  5  - 


as  follows  [21? 


2.12  (D5‘)*  +  §  +  (  ^  -  Eg  a^(a)t)$  =  Q(r) 

2.13  +  (1  -  I  t‘)r^‘  +  (:^t'  »  l)f  =  0 

where  the  apostrophe  represents  differentiation  with  respect 
to  r  (the  radial  distance  from  the  center  of  the  plate), 
t(r)  is  the  thickness  of  the  plate,  D  =  E(t(r))^/12(l  -  v^), 
v(r)  is  the  transverse  deflection  of  the  plate,  crj,(r)  is  the 
membrane  stress  in  the  plate,  q  is  a  uniformly  distributed 
load  acting  transversely  upon  the  plate,  and 

2.14  Q(r)  =  ^ 

2.15  =  aj,(a)t(i') 
and 

2.16  i  = 

Cooi\31nate  system,  gecanetry,  and  variables  are  shown  in  Pig.  3. 

For  the  case  of  the  bending  of  two  joined  cones  the  dif¬ 
ferential  equation  for  each  cone  is  given  by  corresponding  sets 
of  equations  2,3  and  2.4  along  with  definitions  2.5,  2.6, 

2.7,  2.8,  2.9,  2.10  and  2.11. 


>  6  - 


p 


MrofliaiygCT 


IIIo  BOUNDARY  AND  JUNCTION  CONDITIONS 


For  the  cone  with  attached  base -plate  the  boundary  con¬ 
ditions  at  the  free  end  of  the  cone  are 


3.1  M 

at  open  end  of  cone 


and  at  the  center  of  the  base -plate  are 


3.3  $  =  0 

3.4  t  =  1 


at  center  of  base -plate 


where  M  is  moment  per  unit  length  and  e,^  Is  the  strain  of 
*n  ^  ym 

the  middle  surface  of  the  cone  in  the  y  direction  (circum¬ 
ferentially)  . 

In  terms  of  the  variables  of  2.3  and  2.4  the  boundary 
conditions  3.1  and  3.2  may  be  written 


3.5 


En-^  /  Q '  4.  ^  ^ 

12(1  -  v^)  ® 


«  M 


3.6 


h:2h  - 


vh 

B 


)S  + 


P  * 
h'^^S 


at  open  end  of  cone 


At  the  Junction  of  the  cone  and  the  base  plate  the  follow¬ 
ing  conditions  may  be  prescribed 
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3.7 

V 

3.8 

3.9 

«p  = 

ss 

CD 

> 

at  the  junction  of  the 
cone  and  base -plate 

3.10 

"p  = 

®ym 

J 

where  Mp 

and 

tr 

are 

respectively  the 

moment  per  unit  length 

and  the  strain  along  the  circumference  of  the  base -plate. 

In  terms  of  the  variables  of  2.3,  2.4,  2,12  and  2.13  the 
junction  conditions  can  be  written 


3.11  A  =  P  tani,  -  |A_ 


3.12 

12(1  -  v^)  ^ 


3.13  -  f  ==  0 


3 


s 


3.  l4  ^  - 


h(2h'-.  ^)S  +  h^s' 
s 


at  junc¬ 
tion  of 
cone  and 
base¬ 
plate 


Inhere  P  =  Q  =  ^ 


For  the  case  of  the  joined  cones  the  conditions  at  the 
free  ends  are 
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3.15 

3.16 

3.17 

3.18 

and  at 

3.19 

3.20 

3.21 

3.22 

where 


=  M 


1 


M 


^-02 


®yinl  “ 


®ym2  ®2 


the  Joined  ends  are  {[1],  p.  20) 


[(2hh'  -  v|^)S  +  [(2hh*  -  v*!")  S  + 


®1  "  ®2 


-12(1  -r^)*  “'-^1  ^12(1  -  v*")  ® 


I 


(Ptar,(,  .. 


(P  tan,  - 


the  subscrlptB  1  and  2  refer  to  i'’eBpectlve  Joined  cones. 
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IV o  FINITE  DIFFERENCE  SCHEME 

Me  will  approximate  the  differential  equations  for  the 
cone  at  equally  spaced  net  points  s^,  Sg, - ,  Sj  where 


»1+1  -  «!  =  is 

We  will  approximate  the  differential  equations  for  the  base 
plate  by  equally  spaced  net  points  r^,  r^,  - ,  r^  where 


We  will  append  a  subscript  to  each  variable  to  designate  its 
corresponding  approximation  at  a  point.  For  example,  S(s^) 

==  similarly  for  the  remaining  variables. 
Derivatives  will  be  obtained  by  a  Taylor's  series* expansion 
about  surrounding  net  points,  neglecting  terras  0((as)  )  and 


0((&r)‘ 


The  difference  approximations  for  the  cone  equations  are 


4.1 

4.2 
where 

4.3 


^(®i)  *’ii  ^1  ” 

L{sp  + 


p  -  s  V  "J’l-i 

—ss — 


COt'S 


11 


4.4  Pg,  =  [(2  v) 


2Ae 


hj  -i  “  -f-i  “1 

4-  (2s^)  cot^ 

(As)^ 


J 


.2\ 


4.5 


ir-  12(1  -  v'^) 

I'  =  S  Z  COS«fl  SAS  --  s  — — 

i  E  O-  h/  COt^ 


s= 


(Q  )^  COS«J.^ 

^  fO  o  1 
hI^~cot^*^ 


and  the  difference  approximation  for  the  differential  operator 

h  '  U 

4.6  %  -t  (1  +  3s^  “  bT  1 

iL.  •!» 


The  first  and  second  derivatives  are  approximated  as  usual 


by 


4.7 


2  AS 


+  0((as)^) 


and 

4.8 


V. 


i-1 


2V.  +  V .  .  T  o 

1  jkti  +  oatBf) 


{taV 


Substituting  4.7  and  4.8  in  4,6  gives 


L(U^) 


-  {l  + 


3e.j  (h. 


i^“i+l 


h. 


2  AS 


i-1 


-  12 


4.9 


f  2s, 

>  JL 

s 


k]\ 


s. 


'■  (4S)^  24a  J 


1 

1+1 J 


The  difference  approximations  for  the  base -plate  equations 

are 


4.10 


4.11 


II 


J 


«  0 


where 

4.12 


12(1  -  v^) 


Substituting  4.7  and  4.8  In  4.10  and  4.11  gives 


-  13  - 


For  the  case  of  the  cone  with  attached  base -plate  we  first 
investigate  the  difference  approximations  of  the  differential 
equations  at  their  boundaries  and  Jiunction,  l.e.,  where  1  =  0 
and  I,  and  where  j  =  0  and  J. 

For  the  cone  when  1  =  0  the  variables  0  j ,  8^,  0^  j  and 
occur  in  equation  4.1  and  the  variables  S  S^,  and  9^ 

occur  in  equation  4.2.  These  two  relations,  l.e.,  4.1  tuid  4.2 
when  1=0,  will  be  denoted 

4.15  Hite.i.  Bq,  S^)  = 

and 

4.16  S^,  e^)  =  kg 


For  the  base-plate  when  J  -  J  the  variables  $j. 
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fesniirr’-A  "j-i.ijhi.ioi.c.;';;'"  (■ox'por’f:,t;icn 

equation  4.10  and  the  variables  'lij^  tj+i 

occur  in  equation  4.11.  These  two  relations,  i.e.,  4,10  and 
4.11  when  J  -  J,  will  be  denoted 


4.17  ^J4l^  ■'  ^^3 
and 

4.18  ifj, 


In  4.15  and  4.l6  the  variables  and  are  located 
outside  the  junction  of  the  cone  and  base -plate j  similarly  in 
4.17  and  4.18  the  variables  and  tj+i  located  outside 
the  junction  of  the  cone  and  base-plate.  These  four  variables 
can  be  eliminated  by  means  of  the  junction  conditions  3.12  and 
3.14  which  in  difference  form  for  i  =  0,  j  «  J  are 


4.19 


12(1  -  V  ) 


/^J+1  "  ^J-1  .  V  ,  X 

+  a  fj) 


2hv 


E  h^^  8i  -  0  T 

O  +  V 


12(1  -  V  )  2As 


s^  o' 
o 


4.20  j^c:£±l.ll£dL  .  1  .  )  =  h  s  +  h 

S4P  ®  ^  °  °  2i3  ®o  ' 


n.  S-,  -  S  .. 
2  1  ^-1 

’  2as 


Solving  4.15  and  4.l6  for  e_^  and  S__^  and  substituting  in  the 
right  hand  side  of  4.19  and  4,20  respectively  and  solving  4,17 
and  4.l8  for  and  and  substituting  in  the  left-hand  side 

of  4.19  and  4,20  respectively  gives  two  equations  in  the 


variables  e^,  and  S^,  These 

equations  will  be  denoted 


4.21 

“5(80'  V-l)  “  ‘‘5 

4.22 

%(®o’  ®1'  “  ''^6 

VJe  next  eliminate  6^  and  from  4.21  and  4.22  by  means  of  the 
remaining  Junction  conditions  3.11  and  3.13.  In  finite  differ 
ence  form  for  1=0,  J  =  J  these  Junction  conditions  become 


4.23 


ll 


^  P  tailep 


COSw) 


4.24  -  rr  0^ 


Solving  4,23  and  4.24  for  S^  and  9^  respectively  and  substitu¬ 
ting  In  4.22  and  4.21  respectively  gives  two  equations  In  the 
variables  6^,  ^j-1^  V  ^1^  '^J-1'  These  equations  may 

be  denoted 


4.25  63_)  -- 

4.26  ^8^%^  ^1^  ~  ^8 


Now  we  turn  to  the  boundary  conditions.  Setting  1  =  I  in 
equations  4.1  and  4.2  gives  two  equations  which  we  denote  as 


-  16  - 
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4.27  ^9^ ®I+1*  ~  ^9 

4.28  "  ^10 

The  outside  points  and  can  be  eliminated  by  means  of 

the  boundary  conditions  3*5  and  3.6.  In  difference  form  when 
1  =3  I  these  conditions  are 


4.29 


E  h*y»  ©t-jlt  **  6t  I 

_  (Jt±L_lrl  +  ^  G.)  =■-  M 


12(  1  -  v^) 


”  4-1  ^4x0  .  ..  2  ^^I-fl  “  ^I-lx  _  . 

4.30  ^^(2— —  ••  '*■  h^  t - - - —)  ~  e 

2AS  ^  2As 


Solving  4.29  for  and  4.31  ^'or  and  substituting  In 
4.27  and  4.28  respectively  gives  two  equations  which  we  will 
denote  as 


4.31 

H^l(0j, 

®i-i"  4^  ~  4i 

4.32 

42^4^ 

4-1'  ^  42 

Inasmuch  as  the  conditions  3.3  and  3.4  at  the  center  of 
the  plate  are  given  In  terms  of  the  functions  rather  than 
their  derivatives,  the  corresponding  finite  difference  equa¬ 
tions  for  the  base -plate  differential  equations  are  written 
for  J  *=  1 .  The  equations  may  be  denoted 


-  17  - 


4.33 


CvC*  i  i » *  c'v  X  A 1. 1 0  j.  W  O J. ' 


'C.:.  O;  1 


^1'  ^2^  ~  ^13 
4.34  'h*  ^2^  ~  ^l4 


The  boundary  conditions  3.3  and  3.4  In  difference  form  are 

4.35  ® 

4.36  <^0  ^ 

Substituting  4.35  and  4.36  in  4.33  and  4.34  gives  two  equa¬ 
tions  which  we  will  denote  as 


4.37 

*2^  =  *^15 

4.38 

♦2)  “  >^16 

In  4.25,  4.26,  4.31,  4.32,  4.37,  and  4.38  we  have  four 
equations  in  the  four  boundary  variables 
and  two  equations  in  the  two  Junction  variables  §j, 

We  next  obtain  2(l  -  1)  +  2(J  -  2)  equations  at  interior 
points  by  evaluating  4.1  and  4.2  at  1  =  1,  2,  i-i 

and  4.10  and  4.11  at  J  =  2,3,  - ,  J-1. 

When  evaluating  4.1  and  4.2  at  i  =  1  we  bring  in 
and  S^.  In  order  that  all  points  lying  on  the  Junction  satisfy 
the  Junction  conditions,  these  terms  must  be  eliminated  through 
4,24  and  4.23  as  was  done  in  the  case  of  4.21  and  4.22.  Thus 
for  i  --  1  we  can  denote  4.1  and  4.2  as 
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4.39 
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®2^  “  ^17 

4.40  6^^  ~  ^l8 

The  remaining  2(l  -  2)  equations  for  the  cone  may  be  denoted 


4.41 

^li^®l-l^  ®1'  ®i+l' 

Sj)  = 

i  -  2,3 

4.42 

^2i^^l-l»  ®i>  ®i+l' 

=  I'gl 

i  =  2,.  3 

At  the  Interior  points  of  the  base  plate  the  2{j  -  2) 
equations  4.10  and  4.11  may  be  denoted 

4.43  ®3J^^J"1'  ^  ~  j„l 

4.44  ^j+1^  ~  ^^4j  ^  ~  j-1 


We  now  Investigate  the  difference  approximations  for  the 
Joined  cones  at  their  boundaries  and  Junction.  At  the  Junc¬ 
tion,  cone  1  has  its  Ith  net  point  and  cone  2  has  its  0th 
net  point.  Substituting  i  =  0  and  i  =  I  in  4.1  and  4.2  gives 
four  relations  which  we  will  denote 


'(.45  Kj^Ce.i,  e^,  8;^,  Sj,)  =  hj 

t 

4.46  K2(S.j,  s^,  ep  =  hg 


-  19  - 
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4.47 

®I+1*  ~  *^3 

4.48 

^I"  ^1-hi*  Qj)  “  ^4 

In  these  four  equations  we  have  ®I+1'  ^I+l 

lie  outside  the  Junction.  These  four  variables  can  be  elimi¬ 
nated  by  means  of  Junction  conditions  3.19  and  3.21  which  in 
difference  form  are 


4.49 


fhj 


2  AS 


V  h. 


s. 


+  h^^  ( 


^I+l  ■“  ^I-IO 


2^8 


)j 

-*1 


=  Th^  (2 

O  O' 


h,  -  h 


-1 


2  AS 


vh. 


s. 


’)  +  h/  (- 


2  AS 


i)*] 


4.50 


r  ^  ^i^  /®i+i  "  ®i~i  .  V 


E  h  “ 

I®  *0 


h  **  ®-i 


12  (1  -  v'^}  2  AS 


4 


Solving  4.46  and  4.48  for  and  respectively  and  sub¬ 

stituting  in  4.49  gives  one  equation  in  the  variables  S^, 

^1#  Oq,  ^I*  Solving  4.45  and  4.47  foJ*  0_3_ 

and  respectively  and  substituting  in  4,50  gives  one  equa¬ 
tion  in  the  variables  0^,  0^^,  S^,  ^1*  ^I* 
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G* 


1  Techiiology  Coxpoz'-ation 


two  equations  may  be  denoted 


4.51 

KjCS^,  Sj^, 

4.52 

«o'  ^1-1"  hi  “  S 

®o'  hi  **6 


We  next  replace  9^  by  0j  through  Junction  condition  3.20  and 
we  replace  by  Sj  through  Junction  condition  3.22.  Equa¬ 
tions  4.51  and  4.52  then  become 

4.53  K^CS^,,  Sj,  9j)  =  h^ 

4.54  Kq(0^,  ®I-1'  “  ^8 


At  the  free  end  of  cone  1  we  may  denote  equations  4.1 
and  4.2  in  finite  difference  form  as 


4.55  ®1*  ^o^  ~  ^9 

4.56  S^,  S^,  6q)  =  hj^Q 


Boundary  conditions  3.15  and  3.I7  in  finite  difference  form 
are 


4.57 


E  h  3  01  -  9  , 

12(1  -  v'^)  2as  ®o  ° 


4.58 


S,  -  S 


,  ^3  “  ^-1  o  -T  -  -  1 


2as 


)  =  e- 
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Solving  4,57  and  4.58  for  8^-j^  and  respectively  and  sub¬ 
stituting  In  4.55  and  4.56  gives  2  equations  which  we  will 


denote 

4,59 

^11^ ®o^  ®1' 

'•'ll 

4.60 

%2^®o'  ^1'  ®o^ 

=  hj2 

At  the  free  end  of  cone  2  we  have  two  equations  (see 
4.31  and  4.32)  which  we  may  denote 


4.61 

%3^®I'  ®I-1' 

^I^  “  *^13 

4.62 

*^14^®I'  ^I-l^ 

In  4.53i  4.54,  4.59#  4.6o,  4.61^  and  4.62  we  have  four 
equations  In  the  four  boundary  variables  and  of  cone  1 
and  and  S^  of  cone  2  and  we  have  two  equations  In  the  two 
Junction  variables  Bj  and  of  cone  1.  We  next  obtain  2(l  -  l) 
equations  at  the  Interior  points  of  cone  1  and  2(l  -  2)  equa¬ 
tions  at  the  Interior  points  of  cone  2.  When  evaluating  4,1 
and  4,2  at  1  =  1  we  bring  In,  as  mentioned  previously, 
and  S^.  In  the  ease  of  cone  1,  however,  we  already  have  equa¬ 
tions  in  6^  and  S^(l,e.,  4.59  and  4.60)  and  therefore  the 
equations  at  Interior  points  may  be  denoted 
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4.63 

^li^®l-l'  ®i*  ®i+l^ 

sp  = 

4.64 

^2i^^i-l'  ^i'  ^i+l' 

In  the  case  of  cone  2  we  elitolnate  0^  and  through  junction 
conditions  3.20  and  3,22  to  give 


4.65 

^31^®!^  01,  02^  S^) 

■=  hjj. 

4.66 

^4i^®i-l^  ®i'  ®i+l' 

sp  = 

i  ~  2,3,  ^,1-1 

4.67 

^5l(®I>  ^1'  ^2*  ®1^ 

r-i 

li 

4.68 

^6i^^i-l"  ^i'  ^i+1" 

e^)  = 

^6i  ^  “  2,3,  -  -  ,  I-l 

V.  METHOD  OP  SOLUTION  OP  THE  DIPPERENCE  EQUATIONS 

We  can  see  the  form  of  the  difference  equations  for  the 
cone  and  base -plate  equations  if  we  define  the  following 
matrices 
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The  equations  of  the  previous  section  can  thus  be  written  in 
the  form 


The  matrices  A  and  R  can  be  formed  from  equations  4.25,  4.26, 
4.31,  4.32,  4.37,  4.33,  4.39,  4.40,  4.41,  4.42,  4.43,  and 
4.44.  The  matrix  A  is  symmetric  and  its  diagonal  matrices 
are  trl-diagonalj  the  disposition  of  its  elements  is  shown 
on  page  27.  The  equations  corresponding  to  the  elements  of 
each  row  are  given  to  the  left  of  the  matrix.  An  asterisk 
indicates  non-zero  constant  coefficients  and  p  indicates  the 
px’esence  of  the  non-linear  terra  due  to  the  coupling  introduced 


by  2.15.  The  matrices  A^^,  A^g,  Ag^^,  Agg  are  square  with  I 
rovjs  and  columns;  the  matrices  A^^  sind  A^^^^  are  square  with 
J  rows  and  columns.  The  remaining  matrices  are  rectangular 


with  I  by  J  and  J  by  1  rows  and  columns. 


The  solution  of  the  set  of  equations  5.1  could  be  ob¬ 


tained  easily  if  it  were  not  for  the  non-linearity.  Hence 
we  consider  an  iterative  solution  which  assumes  that 


in  2.15  iB  known.  If  we  let  the  superscript  n  designate  the 
nth  iteration,  we  proceed  as  follows:  1)  choose  the  initial 
Cj,(a)  as  2)  solve  the  linear  system  and  calculate  a 

nevt  value  for  a^(a)  from  2.15  using  ij((a),  3)  use  the  new  value 
of  cfj,(a)  to  calculate  a  new  value  of  $(a),  4)  repeat  the 
above  process  until  desired  ©onvergence  is  obtained. 
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We  vjould  expect  this  iterative  scheme  to  converge  be¬ 
cause  is  relatively  small.  If  the  scheme  does  not 

converge,  some  linear  extrapolation  or  Interpotation  between 
succeslve  iterations  may  produce  convergence.  We  leave  this 
question  to  a  future  investigation. 


Finally  we  need  only  give  a  practical  method  of  solving 
the  linear  equations  that  arise  in  the  iterative  scheme. 
f>ince  the  diagonal  matrices  are  trl -diagonal,  we  can  use  the. 
method  of  factoring  the  matrix  into  a  lower  and  upper  matrix. 
This  factorization  requires  the  inversion  of  the 
A^^,  and  A^jj,  ,  The  resulting  equations  can  be  easily  solved 
using  a  digital  computer  [3]. 


In  the  case  of  the  Joined  cones  we  define  the  following 
matrices 


(epi 

Si 

(Sj)o 

^Sl)l 

(^3)2 

Sg  = 

(Sg)^ 

2 
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The  appropriate  equations  can  then  be  written 


The  matrices  B  and  T  can  be  formed  from  equations  4.53, 
4.54,  5.^9,  4.60,  4.6l,  4.62,  4.63,  4.64,  4.65,  4.66,  4.6?, 
and  4.68.  The  matrix  B  Is  symmetric  and  its  diagonal  mati’i- 
ces  are  tri-diagonal;  the  disposition  of  elements  is  shoim  on 
page  28.  The  equations  corresponding  to  the  elements  of 
each  row  are  given  to  the  left  of  the  matrix.  An  asterisk 
indicates  non -zero  constant  coefficients.  Since  all  equations 
are  linear  their  solution  can  be  carried  out  by  a  single 
application  of  the  method  described  previously  for  an  Itera- 
t5-ve  cycle. 
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